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Abstract 

We report new results about the two-time dynamics of an anomalously diffus- 
ing classical particle, as described by the generalized Langevin equation with a 
frequency-dependent noise and the associated friction. The noise is defined by its 
spectral density proportional to uj s ~ x at low frequencies, with < 5 < 1 (subdif- 
fusion) or 1 < 5 < 2 (superdiffusion). Using Laplace analysis, we derive analytic 
expressions in terms of Mittag-Leffler functions for the correlation functions of the 
velocity and of the displacement. While the velocity thermalizes at large times 
(slowly, in contrast to the standard Brownian motion case 5 = 1), the displace- 
ment never attains equilibrium: it ages. We thus show that this feature of normal 
diffusion is shared by a subdiffusive or superdiffusive motion. We provide a closed 
form analytic expression for the fluctuation-dissipation ratio characterizing aging. 
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1. Introduction 

In this paper, we study the two-time dynamics and the aging properties of a 
classical particle diffusing in the presence of coloured (i.e. frequency-dependent) 
noise and friction. The dynamics of the particle is governed by the generalized 
Langevin equation [l]-[2]. We will not be concerned here with the specific micro- 
scopic origin of the noise and of the friction. The knowledge of the noise spectral 
density suffices to make the model precise. The colored noise is conveniently de- 
fined by its spectral density, assumed to be proportional to u/ _1 at low frequencies 
with 0<5<lorl<5<2. The frequency-dependent friction is deduced from 
the noise by means of the second fluctuation-dissipation theorem (FDT). 

This general framework can be used in the classical domain as well as in 
the quantum one, depending on which formulation is chosen for the FDT. The 
one-time dynamics of this dissipative model has been extensively studied, with 
particular emphasis on the quantum case [3] -[6]. In the classical case, the model 
allows for the description of anomalous diffusion, that is, either subdiffusion (for 
< 6 < 1) or superdiffusion (for 1 < 6 < 2). Closed analytic expressions of 
the lowest moments of the velocity and of the displacement of the anomalously 
diffusing particle have been obtained in [7]. 

The two-time dynamics of the diffusing particle has been investigated in the 
case 6 = 1, which corresponds to white noise, non-retarded friction and standard 
Brownian motion [8] -[9]. When studying the two-time correlation functions of 
the dynamical variables attached to the particle, one is faced with interesting new 
physical effects such as aging, i.e. the absence of time-translation invariance. In- 
deed, as first noticed in [8], the two-time correlation function (x(t)x(t')) of the dis- 
placement of the diffusing particle depends explicitly on both times t and t' and not 
only on the difference t — t', even in the limit of large age t' (0 < t' < t). In contrast 
to the velocity, which equilibrates at large times, and does not age, the displace- 
ment is an out-of-equilibrium variable. The equilibrium FDT has to be modified 
in order to relate the displacement response and correlation functions. This can 
be achieved through the introduction of a factor (the fluctuation-dissipation ratio) 
rescaling the temperature [8]. We have demonstrated in [9] that, for a classical 
diffusing particle, the fluctuation-dissipation ratio can be expressed in terms of the 
time-dependent diffusion coefficients D(r) and D(t w ), where r = t — t' denotes the 
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observation time and t w = t' the waiting time. We have subsequently extended 
this result to non integer values of 5 between and 2 [10]. 

In the present paper, we study the two-time dynamics of an anomalously dif- 
fusing particle. Anomalous diffusion, which has been observed in various physical 
systems, is the subject of a growing interest, both experimental and theoretical 
(for a recent review see [11]). Making extensive use of Laplace analysis, we derive 
closed form analytic expressions for the two-time correlation functions of the ve- 
locity and of the displacement. We demonstrate that these two-time quantities are 
expressible in terms of Mittag-Leffler functions [12], which have previously been 
shown to play a central role in the relaxational one-time dynamics of the parti- 
cle [7]. As a result, while the particle velocity thermalizes at large times (slowly, 
except for 5=1), its displacement never attains equilibrium: it ages. We thus 
show that this feature of standard Brownian motion is shared by a subdiffusive 
or superdiffusive particle. We provide closed form analytic expressions in terms 
of Mittag-Leffler functions for the time-dependent diffusion coefficients D(r) and 
D(t w ), as well as for the fluctuation-dissipation ratio characterizing aging. 

2. Generalized Langevin description of anomalous diffusion: a 
Laplace analysis 

2.1. The model and the dynamical variables of interest 

The generalized Langevin equation for a classical particle of mass m diffusing 
in the absence of a deterministic potential writes 

m % + frt- *>(0 dt ' = v = % 

In Eq. (2.1), F(t) is the Langevin force acting on the particle, as modelized by a 
stationary Gaussian random process of zero mean, and ^(t) is a retarded friction 
kernel, which is an even function of t [l]-[2]. As well-known, the coherence of the 
model implies that F(t) and ^(t) are not independent quantities, but instead are 
linked by the second fluctuation-dissipation theorem, that is 

(F(t)F(t')) = kTm 1 (\t-t>\), (2.2) 

where T is the temperature. (The symbol (. . .) denotes the average over the 
realizations of the noise). For white noise, one simply has 

(F(t)F(t')) = 2kTr}8{t-t'), 7 (£) = 2 7 <5(£), (2.3) 
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where 7 is the friction coefficient and rj = rwy the viscosity. The Langevin equation 
is then non-retarded: 

dv „. . . „. 

m — +m^v = F(t). (2.4) 

In the following, we shall be interested with the one-time and two-time prop- 
erties of two dynamical variables attached to the particle, namely its velocity as 
defined by the solution of Eq. (2.1) satisfying the initial condition 

v(t = 0) = v o , (2.5) 

and its displacement, as defined by 

x(t) = [ v{t')dt'. (2.6) 
Jo 

We will primarily be concerned with non-white noise of parameters to be speci- 
fied below and the associated friction coefficient. Let us beforehand describe the 
general lines of our analysis, which heavily relies upon the use of Laplace trans- 
formation. 

2.2. Laplace analysis 

The one-time and the two-time properties of the dynamical variables of inter- 
est can conveniently be studied by means of Laplace analysis, which implies that 
all the quantities of interest are to be considered as causal functions. Defining, 
for any quantity q(t), its Laplace transform q(z) = J °° q(t) e~ zt dt, one gets, by 
applying the Laplace transformation to the generalized Langevin equation (2.1): 

mzv(z) + m A /(z)v(z) = F(z) + mv(t = 0). (2.7) 

2.3. Statistical properties of F(z) 

The generalized Langevin model is fully specified by Eq. (2.7) with the initial 
condition (2.5), together with the given statistical properties of F(z). Since F(t) is 
a Gaussian process, the statistical properties of F(z) are completely characterized 
by the average value (F(z)) and the correlation function (F(z)F(z')}. 

The Langevin force F(t) being a process of zero mean, one has: 



(F(z)) = 0. 



(2.8) 



As for the correlation function (F(z)F(z')), one gets from the fluctuation-dissipa- 
tion relation (2.2): 

/•oo />oo 

(F(z)F(z')) = kTm / e- z,t ' e- zt 1 {\t-t'\)dtdt'. (2.9) 
Jo Jo 

To compute the double integral (2.9), we shall consider separately the two contribu- 
tions (F(z)F(z'))i and {F(z)F(z'))2, coming respectively from the time domains 
t' < t and t < t'. Let us first consider the contribution from the domain t' < t. 
Setting u = t — t', one gets 

/•OO />oo 

(F(z)F(z')) 1 = kTm e- (2+2 ' )t 'cft / / e~ zu 7 (u) du, (2.10) 

that is: 

(F(z)F(z')) 1 = kTm 7(z). (2.11) 

Coming then to the contribution from the domain t <t', one gets 

/•oo />oo 

(F(z)F(z')) 2 = kTm e - {z+z ' }t dt e~ z ' v ^{v)dv, (2.12) 

that is: 

(F(z)F(z')) 2 = kTm i(z'). (2.13) 

z + z 

Adding the contributions (2.11) and (2.13), one finally obtains: 

(F(z)F(z')) = kTm ^ + ^ . (2.14) 

Eq. (2.14), which represents the Laplace domain formulation of the second FDT 
(Eq. (2.2)), will play a central role in the following. 

Actually, relations of the form (2.14) hold quite generally. Indeed, given any 
stationary correlation function of the form 

(<t>(t)<t>(t')) = Af(\t-t'\), (2.15) 

one can show in a similar way that its double Laplace transform writes 

(Hz)Hz'))=aM±M. (2.16) 
z + z 
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2.4. Statistical properties of v(z) 

One deduces from Eqs. (2.7) and (2.5) the expression of v(z): 

fiW = -^r + -^. (2-17) 
As for the average particle velocity, one has, since (F(z)) = 0: 

<*(*)> = JT^Y (2 - 18) 

The double Laplace transform (v(z)v(z')} of the two-time velocity correlation func- 
tion (v(t)v(t')) is given by: 

, HzW)) = + 4 (2 19) 

It is the sum of two contributions, the first one coming from the random force 
correlation function, the second one being linked to the initial condition. Taking 
into account the expression (2.14) of (F(z)F(z')), one can rewrite Eq. (2.19) as 

<OWM) = ^*(4±|m + ^_|:)* w a (i0 , (, 20 , 

with K(z) as defined by 

K{z) = 1 —- . (2.21) 

z + -f(z) 

Formula (2.20) displays clearly the fact that the correlation function (v(z)v(z')) is 
the sum of two contributions of a different character as far as stationarity prop- 
erties are concerned. The first one, ~ ^~^~t^^ ? of the general form (2.16), 
corresponds in the time domain to a stationary random process. The second one, 
(vq — ^) K{z) K(z'), corresponds to a function depending separately on t and t' 
(and not only on the difference t — t'). This term, if non-zero, will yield an aging 
contribution to the two-time velocity correlation function (v(t)v(t')). 

Averaging out Eq. (2.20) over an equilibrium ensemble of initial velocities, 
the aging contribution to (v(t)v(t')) disappears. We shall come back to this point 
later. One is then left with: 

m^))= k - k(z)+k{Z,) . (2.22) 
x ' " m z + z' v ' 
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2.5. Statistical properties of x(z) 

One deduces from Eq. (2.6) the expression of x{z): 



x(z) = -v(z). (2.23) 

z 

As for the average particle displacement, one has: 

<*(*)> = 7 T^M- (2-24) 
2^ + 7(2j 

We assume that the particle velocity is at thermal equilibrium and does not age. 
The double Laplace transform (x(z)x(z')) of the two-time displacement correlation 
function (x(t)x(t')) is given by: 

i z x z' = — ^ 2.25 

x w v ;/ m z z' z + z' v ' 

Due to the presence of the factors 1/z and 1/z' in the r.h.s. of Eq. (2.25), 
the correlation function (x(z)x(z')) is not of the form (2.16). Consequently, in 
the time domain, the two-time displacement correlation function {x(t)x(t')) will 
depend separately on t and t'. The particle displacement is thus an aging variable. 

Let us now be more specific and introduce the parameters of the so-called 
non-Ohmic noise and friction allowing for the description of the dynamics of a 
subdiffusive or superdiffusive particle. 

3. Noise and friction of the non-Ohmic Langevin model 

We consider a coloured noise of spectral density 

<|FM|'> = 2*T W (H)*- 1 /«(H). (3.1) 

The small- 1 a; | behaviour of (|-F(u;)| 2 ) is a power-law characterized by the exponent 
5 — 1. The function f c (\uj\/uj c ) is a high-frequency cut-off function of typical width 
co> c , and u> <^ u c denotes a reference frequency allowing for the coupling constant 
rjs = rwys to have the dimension of a viscosity for any 5 [6]. For 6 = 1, the 
noise spectral density is a constant (white-noise), at least in the frequency range 
\u\ <^ uj c . Then the Langevin force is delta-correlated and the Langevin equation 
is non-retarded. The white-noise case corresponds to Ohmic friction. The cases 
< 6 < 1 and 5 > 1 are known respectively as the sub-Ohmic and super-Ohmic 
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models. Here we will assume that < S < 2, for reasons to be developed below 
[3]-[6]. 

To begin with, we have to compute the correlation function (F(z)F(z')). We 
first deduce from Eq. (3.1) the expression of the two-time correlation function 
(F(t)F(t')) by using the Wiener-Khintchine theorem: 

(F(t)F(t')) = 2kT V r ^(M) 5-1 /c (M) cosc^-O- (3-2) 

Then, one deduces from Eq. (3.2) the expression of the double Laplace transform 
{F(z)F(z')): 



z + z' J tv \uj/ W c / 



+ 



Z 2 + UJ 2 z' 2 + UJ 2 



(3.3) 

This expression is actually of the form (2.14), with 



= 15 



r-w^wji) * (3 . 4) 

J TV \UJ/ Vcu c / Z Z +UJ Z 



Since < 5 < 2, one can safely make the limit u c — > oo in formulas (3.3) and 
(3.4), and write 

du/un 5 - 1 z 
n \u)J z 2 + u) 2 

and 



7(*) = 7*/ -(-) (3-5) 



(F(z)F(z')) = 2kT m - J— 1 f -^-^ + -^U^ • (3.6) 

The integrations over a; in formulas (3.5) and (3.6) once carried out, one gets 

^)=^(£f _1 ^' (3 - 7) 

W/ sin ^ 

and 

At this stage, following [6], it is convenient to introduce the 5-dependent 
frequency us as defined by 
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Using this notation, one has 



J 5 

and 



7(z)=wJ-V- 1 , (3.10) 



(F(z)F(z')) = kTmu 2 s ~ 5 ~ ^ + * ~ ) . (3.11) 

Eqs. (3.10) and (3.11) together with the Langevin equation (2.7) and the initial 
condition (2.5) fully specify the model in the Laplace domain. 

4. Particle velocity: one-time and two-time properties 

4.1. The average velocity 

By inverting the Laplace transformation, one gets from Eq. (2.18): 

i pc+ioo i 

WW^L. 7TW) e " dz - (41) 

The constant c in the definition of the integration contour in Eq. (4.1) is real and 
chosen such that all the singularities of K(z) = l/[z + 7(2)] are lying to the left 
of the integration path. Using the expression (3.10) of 7(2), one is left with 

W))=v«— (4.2) 

Ain J c — loo Z -+- UJg Z 

i.e., as showed in [6]-[7]: 

(v(t))=v E 2 - 5 [-(u s t) 2 - 5 ]. (4.3) 

Eq. (4.3) displays the fact that the average velocity relaxes towards zero, this 
decay being described by the Mittag-Leffler function 1 E a (x) [12] with a = 2 — 5 
and x = -(uj s t) 2 - s [6]-[7]. 



1 The Mittag-Leffler function is defined by the series expansion 



««M = E rjJVi). - > 0. («) 

n=0 v 7 

where V is the Euler Gamma function. The Mittag-Leffler function E a (x) reduces 

to the exponential e x when a = 1. The asymptotic behaviour at large a; of the 
Mittag-Leffler function E a (x) is as follows: 

Ea(x) ^ -- r , X » 1. (4.5) 

x 1 (1 — a) 
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At large times (i.e. oj$t 3> 1), the average particle velocity decreases according 
to a power-law of time: 

{v(t))~volj0£^, 0<5<2, 8 + 1. (4.6) 

In the Ohmic case 8=1, one has ujs=i = 7- Eq. (4.3) is then nothing but the 
standard exponential decay of the Brownian particle average velocity: 

(v(t)) = v e-^. (4.7) 



4.2. The velocity correlation function 

Eq. (2.20) yields in the time domain 

kT 

( v (t)v(t')) = —E 2 - S [-Mt-t'\ 2 - 5 )] 
m 

+ (v 2 - E 2 _ 5 [-(ivst) 2 - 5 } x E 2 _ 5 [-(ivst') 2 - 5 } ■ (4-8) 

As expected, the two-time velocity correlation function (v(t)v(t')) is the sum of a 
stationary part and of an aging one. In the Ohmic case 8 = 1, Eq. (4.8) reduces 
to: 

(v(t)v(t')) = — e-Tl*-*'l + (v 2 - —) e~^ t+t '\ (4.9) 
m \ m J 

As already noted, the aging part of {v(t)v(t')) vanishes if one carries out an 
average of Eq. (4.8) over an equilibrium ensemble of initial velocities. Equivalently, 
this can also be achieved if the initial time, instead of being taken equal to 0, is 
taken equal to U, and subsequently rejected to — oo. The aging part of (v(t)v(t')}, 
which is then equal to the product E 2 s(-[uJs(t-ti)] 2 ~ s ) x E 2 s(-[uJs(t' —U)] 2 ~ S ) , 
indeed vanishes in the limit £j — > — oo. In this limit, any initial fluctuation of the 
average velocity has decayed to zero at finite time. The particle velocity is then at 
equilibrium. Correspondingly, the two-time velocity correlation function reduces 
to its stationary part: 

kT 

(v(t)v(t')) = —E 2 _ s [-(oj s \t-t'\) 2 - s )}. (4.10) 

Interestingly enough, the stationary part of (v(t)v(t')) is proportional to the 
Mittag-Leffler function E 2 _ 5 [-{u 5 \t - t'l) 2 " 5 )] (Eq. (4.10)), while the average 



11 



velocity (v(t)) resulting from a given initial fluctuation evolves proportionally to 
E2-5 [— (uJst) 2 ~ s ] (Eq. (4.3)). This fact constitutes the generalization to the non- 
Ohmic case (5 7^ 1) of the regression theorem according to which the fluctuations 
decay in time following the same law as the average value, as it is well-known for 
the exponential decay of fluctuations which holds in standard Brownian motion 
(5 = 1). 

In the particular case t = t', one gets from Eq. (4.8) 

<^)> = ^+(-o 2 -^) {E 2 - 5 [-^stf- 5 ]}\ (4.11) 

a result which describes the slow relaxation of the second moment of the velocity 
towards its equilibrium value [7]. 

5. Aging of the particle displacement 

From now on, we assume that the two-time velocity correlation function re- 
duces to its stationary part (4.10). The particle velocity is thermalized. Then, the 
two-time displacement correlation function can directly be obtained by Laplace 
inversion of the correlation function (x(z)x(z')) as given by Eq. (2.25). One thus 
gets, for < t' < t: 
kT r 

(x(i)x(t')) = — \t 2 E 2 _ s ,s[-(^t) 2 - S } +t' 2 E 2 _s4-(^t') 2 - S ] 

~(t-t'fE 2 ^{-[u 5 {t-t')] 2 - s )), (5.1) 

where we have used the generalized Mittag-Leffler function 2 E a ^(x) [12]. Another 
way of deriving the result (5.1) is to compute the double integral 

(x(t)x(t'))= I [ (vit^vih)) dhdh (5.3) 
Jo Jo 

the correlation function (v(ti)v(t2)) being taken from Eq. (4.10). One has to use 
the integration formulas 

/ E 2 _ s [-(u 5 t 1 f- 5 ]dt 1 =tE 2 ^ 2 [-(u s tf- 5 ] (5.4) 
^0 



The generalized Mittag-Leffler function is defined by the series expansion 



Tl 
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and 



[ dh I' E 2 - S [-(u s t 2 ) 2 - s ]dt 2 =t 2 E 2 - S , 3 [-(ust) 2 - 5 }. (5.5) 
Jo Jo 

Eq. (5.1) displays the fact that the particle displacement is not a variable at 
equilibrium. When 8 = 1, we recover the result we previously obtained in [9]: 

kT f 1 + e- 7(t_t,) - e _7t - e -7 *' \ 

(x(t)x(t')) = — \ 2t'- — . (5.6) 

V \ 7 



5.1. The equal time correlation function and the time-dependent 
diffusion coefficient 

In the particular case t = t', one gets from Eq. (5.1): 

kT 

(x 2 (t)} = 2 — t 2 £ 2 _ 5 , 3 [-(ojst) 2 - 6 } . (5.7) 
m 

The time-dependent diffusion coefficient, as defined by [9] 

m = \jt {x2{t)) ' (5 - 8) 

is: 

kT 

D(i) = —tE 2 _ 5 , 2 [-{uj s t) 2 - 5 ]. (5.9) 

In the Ohmic case 8=1, Eqs. (5.7) and (5.9) reduce to the standard formulas 

kT 

(x 2 (t)) = 2 — t (5.10) 
my 

and 

kT 

D(t) = —(l-e-^). (5.11) 

my 

At large times (i.e. uj§t 3> 1), Eqs. (5.7) and (5.9) yield 



kT 1 (u 5 t) s 
~m T(6 + l) 

and 

kT 1 (cost) 6 - 1 



(*\t))^-.T^ (5-12) 



m 005 1 (o) 

Eqs. (5.12) and (5.13) show that, for oj$t 3> 1, the particle motion is subdiffusive 
for < 8 < 1 and superdiffusive for 1 < 8 < 2. In view of this fact, the expression 
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(5.1) of {x(t)x(t')) shows that the displacement of a subdiffusive or of a superdif- 
fusive particle is an aging variable, as it is the case for a normally diffusing particle 
[8]-[9]. 

5.2. The fluctuation-dissipation ratio 

Concerning the particle displacement, which is not a variable at equilibrium, 
a modified FDT can be written as 

Xxx (t,t') = m - )} , (5.14) 

where Xxx(t,t') is the displacement response function [8]. We have previously 
demonstrated that, for a diffusing particle, the fluctuation-dissipation ratio X(t, t') 
can be expressed in terms of the time-dependent diffusion coefficients D(t) and 
D(t w ), where r = t — t' denotes the observation time and t w = t' the waiting time 
[9]-[10]: 

Using the above found expression (5.9) of the time-dependent diffusion coefficient, 
one gets: 

When 6 = 1, we recover the result we obtained in [9]: 

X(r, t w ) = )' e T 7T . (5.17) 

At large observation and waiting times (i.e. ujst ^> 1, ujgt w ^> 1), one can 
use in Eq. (5.16) the asymptotic expressions of D(r) and D{t w ) as given by 
Eq. (5.13). Eq. (5.16) then displays the fact that, in a subohmic or superohmic 
model of exponent 5, the large times aging regime is self-similar, as pictured by 
the fluctuation-dissipation ratio [10]: 

X(r,t w )~ j!— j_ (5.18) 

1 + (t) 

Interestingly enough, X(r, t w ) is then a function of t w /r, solely parametrized by 
5. For 6 = 1, one retrieves the result X = 1/2 [8] -[9]. For any other value of 5, X 
is an algebraic function of t w /r. The limits r — > oo and t w — > oo do not commute. 
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6. Conclusion 

We have studied the two-time dynamics of a classical anomalously diffusing 
particle. The particle motion is governed by the generalized Langevin equation, 
with a noise of spectral density proportional to uj 5 ~ 1 at low frequencies (0 < 5 < 1 
or 1 < 5 < 2). The Langevin force is modelized by a stationary Gaussian random 
process. The velocity and the displacement of the particle are themselves Gaus- 
sian processes, fully characterized by their averages and their two-time correlation 
functions. 

Using a double Laplace transformation technique, we have been able to derive 
closed form analytic expressions for the correlation functions of the velocity and 
of the displacement. Both can be expressed in terms of Mittag-Lefffer functions, 
which are thus demonstrated to play a central role, not only in the one-time 
dynamics as shown in [7], but also in the two-time dynamics of the anomalously 
diffusing particle. While the velocity thermalizes at large times (slowly, except 
for the standard Brownian motion case 5 = 1), the displacement never attains 
equilibrium and ages. 

The thermal equilibrium velocity fluctuations decay in time following the same 
Mittag-Lefffer law as the average velocity resulting from a given initial fluctuation. 
This result constitutes the generalization to the non-Ohmic case 5 ^ 1 of the 
regression theorem well-known for the exponential decay of fluctuations which 
holds in standard Brownian motion. 

As for the particle displacement, we have obtained for the fluctuation-dissipa- 
tion ratio characterizing aging an expression in terms of Mittag-Lefffer functions, 
valid for any values of the observation and waiting times. For < 6 < 1 or 
1 < 5 < 2, the large times aging regime is described by a self-similar function 
of t w /r. This is in marked contrast to the standard Brownian motion case, in 
which the large times aging regime is described by the value X = 1/2 of the 
fluctuation-dissipation ratio. 
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